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ABSTRACT 


Generation of parametric x-radiation (PXR) may be 
described as the Bragg scattering of virtual photons to 
produce real x-rays which satisfy the Bragg condition nd = 
2dsin@,, where o., is the angle between the electron beam and 
the crystal plane. Enhanced higher order parametric x- 
radiation from the <220> and <111> planes of silicon 
crystals of varying thicknesses were observed. Production 
of PXR of the order n=1 for both planes of a 20 Um thick 
crystal and of the orders n=1, and n=2 of the <220> and the 
n=1, n=3, and n=4 of the <111> planes of a 44 Um and 320 um 
crystal were observed. Exploiting the formation. and 
attenuation lengths oof silicon crystals of various 
thicknesses, higher order x-radiation production is enhanced 
relative to the lower energy first order x-ray. Photons of 


4.5 to 21 keV have been observed. 
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I. INTRODUCTION 


The need for a monochromatic and intense x-radiation 
source has led to the study of the interaction of 
relativistic charged particles in various media. Parametric 
x-radiation (PXR) is a possible source.[Ref. 1] PXR is 
formed by a relativistic charged particle, an electron in 
this case, interacting with the planes of a crystalline 
structure at or near the Bragg condition. Ter-Mikaelian 
first theorized this phenomenon for the number of x-rays 
produced ina thin crystal[Ref. 2]. In 1985, Baryshevsky et 
al conducted the first experiment at Tomsk to observe PXR. 
900 MeV electrons bombarded a diamond crystal at a Bragg 
angle of 45 degrees producing parametric x-rays in the 5 to 
25 keV range.[Ref. 3] 

The Naval Postgraduate School (NPS) experiment was the 
first measurement of PXR in the United States. The PXR 
source crystal for the experiment described in this thesis 
was silicon. 

The measurements in this thesis demonstrate how the 
enhancement of higher order PXR due to the x-ray formation 
and attenuation in silicon were affected by the crystal 
thickness. This thesis continues the investigation started 


in Reference 4. 


—— 


II. THEORETICAL BACKGROUND 


Parametric x-radiation is one of several methods for 
producing x-rays when a charged particle beam interacts with 
a crystal. Other production methods include transition, 
diffraction, and channeling radiation. Comparison of the 
various x-ray sources has shown PXR gives the highest 
spectral intensity.[Ref. 5] 

When a charged particle beam satisfies the Bragg 
condition with the crystal lattice, the virtual photons 
associated with the relativistic particles are diffracted. 
PXR is produced when the periodic electric susceptibility 
associated with the crystal lattice structure interacts with 
the field of the particles.[{Ref. 2] 

The periodic arrangement of atoms in a crystal leads to 
the diffraction of x-ray photons with the wave vector k 
which satisfies the Bragg condition 


2 
k* = (kg +t), Ik| = |Ko 





' (1) 
where T is the reciprocal lattice vector and Kp is the wave 
vector of the incoming virtual photons. This process is 
illustrated in Figure 1. In the x-ray range, the refractive 
index n(K,@) becomes greater than unity and the Vavilov- 


Cherenkov radiation requirement, 





k : 
pe coy) eee ) 
@ 


is fulfilled even though the dielectric constant is less 
than one. PXR originates when both the Bragg and Vavilov- 
Cherenkov conditions are satisfied by the wave vector k and 
the frequency @ of the photons generated by the charged 
particles moving through the lattice structure of a crystal 
with uniform velocity.[(Ref. 3] In contrast to the crystal 
structure, homogeneous media require that the refractive 


index n be less than one for the production of x-rays, thus 





Figure 1. Diffraction of virtual photons associated 
with a relativistic charged particle by the planes 
of a crystal [Ref. 2]. 


demonstrating that under the Bragg diffraction condition the 


Vavilov-Cherenkov radiation is PXR. 


Kinematic theory of PXR is applicable for real crystals, 
which consist of thin mosaic blocks turned relative to each 
other at the angle 6 = m/E. Here m is the mass of the 
particle in the beam and E is the energy of the particle in 
the beam. Ter-Mikaelian formulated the expression for the 
number of photons emitted by a charged particle in a thin 
crystal; treating the interaction as a perturbation with the 
crystal requirement 

kL|n-1| < 1, (3) 
where L is the thickness of the crystal. If the refraction 
of photons and the multiple scattering of electrons in the 
crystal are accounted for, then the results, for a mosaic 
crystal with 5 > m/E and with PXR from different blocks 
being incoherent, remain correct. The spectral and angular 


distributions of emitted photons are formulated by 








ON _ e? 2 
On, dw _ on Wald, | 


Cen - Mag), (4) 


q = < ~ Vw? ~ x? + 12 - “[Re(go) - 48], (5) 


~ EVER (6) 


The Fourier components |g;| of the dielectric susceptibility 


corresponding to the reciprocal lattice vector and the mean 
value g, are directly related by the coherent scattering 
amplitude of the photons. 8, is the angle of multiple 
scattering, E, = 21 MeV, and Lp is the radiation length. 
Equation (4) is written in a coordinate system with the z- 
axis directed along the particle velocity v; ky) = On); and 
La = [wIm(g,)]~?. (Ref. 5] 

The characteristic of producing quasimonochromatic 
photons at large angles relative to the velocity vector, v, 
distinguishes PXR from other radiation mechanisms. The 
particle energy does not determine the emitted photons wave 


vectors, k, nor frequencies The photons are defined by 


B° 
the reciprocal lattice vector, T, and the Bragg angle, 05: 


between v and the crystallographic planes associated with T. 


The PXR maxima are defined by the crystal lattice structure 


amplitudes corresponding to tT. The values of the emitted 
wave vectors, k, and frequencies @, are 
Opv k-v 
i ee —— = cos20,, (7) 
Vv Op 


T NIC 


eae os fs a n=1,2,3... 8 
2sin@gp dsin@,’ = (8) 


r2 
Wp(n) = Ata = 


where d is the distance between the crystal planes 


corresponding to the vector T.[Ref. 3] 


The intensity of PXR is proportional to |g,|@ and 


therefore it is greater in crystals with higher packing 
factors, such as a diamond lattice. Rotation of a crystal 
by some A@ results in the emitted wave vector, k, rotating 
2A8 in the same direction while the frequency varies as a 
result of equation (8). The angular distribution of emitted 


photons can be derived by integrating equation (4) over QW, 


resulting in 





g, (wp(n))|* [0% cos” 26 + 09] 


ae - Ya 
a 1 — ; 
06,06, pipe Bal : | sin’ 6p [ 6% A 0g a 631] e 
(k 7 Ky 5 m? 
Oxy = ic Spn= Szesi= Igo. (10) 


where m is the rest mass of the electron. 


The frequency distribution can be determined by 


integrating equation (4) over ny. resulting in 


ON 
ae = Ny (u); 


1 + u*(1 + cos” 26p) _ 


J z(u) =  (aeuye ' (11) 


sin@p (@ - wp) 
cos6g Oph 





T 
Meee 0" u= 


The number of # photons’ detected in an-— angular 


distribution of oF about 0, is given by 


3 
Np = TN o( 1 <i cos” x). "Gta? = 


(12) 
06 + 621, 66, | 
= Nz\1 + cos’ 20 | = , 
| an” B+ Or 
where 
8p 
Li Ss 13 


The detector angular size is a factor in determining the 
value of Nye even for 8, >> eane This factor does not occur 
in channeling or in bremsstrahlung radiation; however in 
PXR, the intensity slowly decreases as 91 away from the 
Bragg condition angle. 


The angular distribution as defined by the parameter 


Bon (®) increases for the electron energy 


(14) 





thus setting the limiting threshold energy for PXR. The 


amount of photons decreases rapidly for E < Eo The 


DE: 
number detected becomes 


4 4 
6 E 
Np = : No(1 te cos’ 205) 3( 22 = ea . (1.5) 


This phenomena illustrates the problem that--even if all 
other parameters are met--if the beam energy is not 


sufficient PXR may not be observed.[Ref. 5] 


Reference 4 was used extensively for the presentation of 


the derivation given above. 


IItl. PXR EXPERIMENT 


The dependence of parametric-x-ray generation on crystal 
thickness was studied using silicon crystals with 
thicknesses of 20 Um, 44 Um, and 320 Um. Two specific 
properties were examined; first, how was the PXR intensity 
effected by the crystal thickness, second, how did the PXR 


linewidth change as a function of x-ray photon energy. 


A. EXPERIMENT DESCRIPTION 

The NPS Linear Accelerator (LINAC) was used to conduct 
this experiment. Since the LINAC is pulsed at 60 Hz with a 
pulsewidth of one microsecond, the low duty factor required 
the measurements be made with only dark current, estimated 
to be much lower than a nanoampere average current. The 
electron beam energy was 85 MeV. 

The experimental set-up is pictured in Figure 2. live 
consisted of a scattering chamber with a target support 
ladder (holding the silicon crystal) in the vacuum chamber 
center. An alignment/reference laser was located outside 
the vacuum chamber 90 degrees counter-clockwise from the 
incoming particle beam inlet. Attached to a chamber port 45 
degrees clockwise from the beam exit port was a brass pipe 


fitted with a one mil Kapton window (1" in diameter) at the 


end immediately before the Si(Li) x-ray detector. 


The 


detector is an EG&G ORTEC system, model number 7113-16250. 





Figure 2. End station PXR experimental set-up. 
1.Incoming electron beam. 2.Target ladder with silicon 
crystal and calibration foils. 3.Scattering chamber. 
4.  Reference/Alignment laser. 5. Vacuum coupled 
flange. 6. Si(Li) x-ray detector. 


Located on the target ladder was a silicon crystal. 
Spacing ,d, of the <111> planes is 3.135 angstroms while 
spacing of the <220> planes is 1.920 angstroms. Silicon 
a density of 2.33 grams per cubic centimeter. For 


experiment two separate crystals were used. The first 


The 


the 


has 


the 


was 


etched to a thickness of 20 tim in a one square centimeter 
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Square in the center of the crystal. The second crystal was 
etched to 44 um in a one square centimeter Square at the 
center. The thickness of the area outside the center of the 
latter crystal was 320 Um. The tolerance in the thickness 
measurements were unknown, and the thickness uncertainties 
were assumed to be negligible for the main conclusions of 
this study. 

8. is the angle between the alignment laser and the back 
surface of the crystal and is used for positioning the 
appropriate planes of the crystal at the Bragg condition 
angle, 0, - The laser defined the zero or home position. 
The incoming beam entered the scattering chamber through the 
270° port. The Bragg angle, On: was referenced to the 
electron beam. as indicated by Figure 3. The Bragg angle 


for this experiment was 22.5°, with the detector located at 


an angle of 9, = 20, - 


B. ENERGY CALIBRATION PROCEDURE 

To calibrate the parametric x-ray energy, a calibration 
target consisting of stacked Cu and Ti foils that was 
overlayed on the backside of a phosphorescent’ screen. 
Figure 4 is a drawing of this device. 

PXR must be detected during the on-time of the electron 
beam pulse. Consequently, the detector and its associated 
electronics exist in an environment which is affected by the 


klystron RF noise and ground loops. Measurements of x-rays 
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Alignment Avs 
LASER Si(L1} 
Detee- or 





Figure 3. This picture depicts the relationship 
between the incoming charged particle beam and the 


crystal planes and the angle, 90,., between the 
alignment laser and the crystal's back surface. 
The detector angle, 8, is twice the Bragg angle, 


0, - 





(D) | 
Figure 4. The drawings indicate how the calibration 
foils were placed on the phosphorescent screen. (4@). 


is the backside of the screen. (b) is the left edge 
of the screen. 
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during the beam pulse are affected by the noise environment 
and will appear at apparent energies which do not correspond 
to the calibration from standard sources taken when the beam 
was off. Consequently a technique had to be devised to 
energy calibrate the spectrum during the machine on-time. 
During the course of this experiment a method was developed 
to calibrate the detector using x-ray floresence. 

The procedure required that the Cu and Ti foils be 
placed in the path of the incoming electron beam. De- 
excitation K-edge x-rays following the interactions of the 
electrons in the beam with the atoms in the foils were 
observed. The x-rays that reached the Si(Li) detector were 
counted at a rate of one photon per machine pulse. (The 
accelerator pulse at a rate of 60 Hz with a pulse duration 
of one microsecond. ) 

With a spectrum of titanium and copper Kg peaks--signal 
to noise ratio on the order of 10 to 1--the resolution of 
the detector could be approximated. The resolution of the 
detector could then be applied to the PXR spectra. 

The Kg x-ray peak of titanium should be observed at 4.5 
keV while that of copper should be seen at 8.05 keV [Ref. 
ile Figure 5 is a plot of the calibration spectrum taken 
July 22, 1991. In it, the signal to noise ratio for both 


the titanium and the copper peaks is approximately 10:1. 


Additionally, the Kg and the Kg lines were resolved at 8.05 
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keV and 8.94 keV respectively. On September 05, 1991, 
another calibration spectrum was made (Fig. 6). Here the 
Signal to noise ratio of titanium was approximately 3:1, 
while that of the copper was closer to 2:11. During the PXR 
experiment, only the Ti and Cu calibration lines were 
obtained. In order to obtain a calibration line for all the 
peaks observed, a subsequent measurement of the fluorescent 
x-ray from titanium, copper, yttrium and indium was made. 
These are shown in Figures 7 thru 10. Figure 11 is a 
composite spectrum created from the data of Figures 7 to 10. 
This spectrum allows for calibrations of detector spectrum 
energy and width for the full range of x-rays from 5 to 30 
keV. A follow-on thesis will develop a sandwich foil which 
will allow simultaneous measurements of all fluorescent 


lines in this energy range. The Kg x-ray peaks for yttrium 


and indium are 14.920keV and 24.104keV respectively while 
the KB peaks are at 16.87 keV and 27.57 kev [Ref. 6]. With 
the x-ray energy peaks from the last set of de-excitation 
data, the detector resolution could be determined. 

To determine the detector resolution, the calibration 
peaks were assumed to be Gaussian in shape. Next the 


natural linewidths of the foils had to be determined. Using 


IThe resolution of the titanium peaks Is approximately the same when uncertainties 
are neglected. 
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1991, using Ti and Cu foils in dark current. 
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Calibration spectrum taken on September 
with Ti and Cu foils in dark current. 
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1992, of Y, foil in dark current. 
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Figure 10. Calibration spectrum taken April 30, 


Le92, 


of In foil in dark current. 
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Figure 11. Composite calibration spectrum taken 


from April 30, 1992, Ti, Cu, Y, and In foils data 
(Figures 7 to 10). 
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Reference 7, the natural linewidths varied in size from 10 
ev to 45 ev. The natural linewidths were also assumed to 
be Gaussian distributions. As both peaks were assumed to be 
Gaussian in nature, the observed linewidth is related to the 
detector resolution and natural linewidth by 

Téps= TGett fw (16) 
where Iopc6 is the measured resolution, Iget is the detector 
resolution and [yy is the natural linewidth. 


Since the uncertainty in I}, 1s assumed to. be 


negligible, the uncertainty in Iget is the same as that for 
lobs: 

Table 1 lists the values ef the Kg x-ray energy peaks, 
the foil natural linewidths and the observed linewidths? for 
the July 22, and April 30, spectrums. Table 2 contains the 
results of using equation (16) to calculate the detector 
linewidth. Also in Table 2 is a listing of the theoretical 


detector resolution as determined from 


— 


Taet= V(ln)* + (2. 35VEFE)*, a! 
where IT, = 110 ev, € =3.8 ev, F, the Fano Factor is 0.125, 
and E is the x-ray energy [{Ref. 10]. Note the manufacturers 
theoretical values do not agree with measurements. Thus the 
value of the detector resolution to be used will be that of 
the observed calibration spectrum linewidths. 
2Three independent methods were used to determine these values. They included 


using Personal Computer Analyzer (PCA) software [Ref. 8], spreadsheet software [Ref. 9]. 
and measuring the linewidths by hand. 
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Figure 12 is a plot of the detector resolution versus 


peak energy for the calibrations done on July 22, and April 










| dan (Ka) | 24.10 .0437 ow. 7007 
fim exp fovisr | + | » | * |e 


Table 1. List of the calibration foils, the respective x- 
ray energy peak values. Included in the list are the 
natural linewidths of the foils and the observed linewidth 
values determined from _ PHA, spreadsheet, and hand 
measurements. 


a0. Of particular interest is the differrence in the two 
titanium widths and the two copper widths. Due to this 


difference and the lack of reference energy peaks in the 15 


3[Ref. 8]. 

4Hand Measurements. 

S[Ref. QI. 

6The linewidth of Y wasn't listed in Reference 7,. The value was approximated by using 
the values of Sr and Zr which was .035 for both. 

’The linewidth of In was not listed in Reference 7, but was approximated from Pd and 
Ag which are both at .043. 
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keV and 25 keV region in the July 22, data, a least squares 


fit was applied to the April 30, spectrum to determine the 


Date Element | Energy | Tae? Vaet? Taett? | Taet 
(kev) (keV) (keV) (keV) (kev) 
amet, fasion {a 






Table 2. List of the Sceeceor a as aeCernnnieS by 
PCA, spreadsheet, and hand measurements. The last column 
represents the calculated resolution from Reference 10. 


slope. The resolution as a function of photon energy was 
assumed to be linear. The linear regression was done using 
an HP15C hand held calculator. The slope was m = 0.0068 
with an intercept of b = 0.543 keV and a correlation r = 
0.978. Next, the value of the slope was used to plot a line 
to best approximate a straight line fit of the July 22, 
data. Figure 13 is a plot of the 22 July, data and the 
manufacturer's theoretical values from Table 2. Note that 
the manufacturer's data plots as nearly a straight line. 


The slope is m = 0.006, the intercept is b = 0.131 keV, and 


8iRef. 8]. 
Hand Measurement 
10iRef. Q]. 
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Figure 12. The plot depicts detector resolution 
versus foil peak energy for the calibration foils. 
x represents the July data while + represents the 
April data. The linear regression is depicted by 
a solid line and the transposed line used for 
translation of the PXR data is represented by the 
dashed line. 
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Figure 13. The plot depicts detector resolution 
versus foil peak energy for the calibration foils. 
x represents the July data while + represents the 
manufacturer's data. The linear regression is 
depicted by a solid line and a translated line 
used used to approximate a straght line fit to the 
data is represented by the dashed line. 
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the correlation is r = 0.9967. Consequently, the assumption 
that the resolution increases linearly aS x-ray energy 
increases is valid. Thus, the transposed line in Figure 12 
would then define the detector resolution for the PXR 


spectrums. 


C. PXR DATA COLLECTION 

With the initial calibration spectrum completed, the 
target ladder was lowered to a position such that the center 
of the 20 [Lm silicon crystal was aligned for the incoming 
electron beam. The ladder was then rotated to 6, = 67.5°!1!, 
This angle corresponded to a Bragg angle of 22.5° for the 
<220> planes of the crystal. A single peak corresponding to 
the n=1 peak at an energy of 8.63keV was observed as shown 
in Figure 14. Following this run, the ladder was rotated to 
the reference point and then rotated to 8, = 102.75°. In 
this position, the <111> planes of the crystal were aligned 
at the 22.5° Bragg angle. This result is plotted in Figure 


15. The n=1 peak was observed at 5.19 keV. 


11 The silicon crystals were provided by Adelphi Technology Inc., and were etched so 
the <220> plane was perpendicular to the incoming electron beam when 86, = 0: 
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Figure 14. The plot is the n=1 peak of the 20 [im 
Si <220> orientation using an 85 MeV dark current 
electron beam. The peak is located at 8.63 kev. 
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Figure 15. The n=1 peak of the 20 Um Si <111> 
Orientation produced by an 85 Mev dark current 
electron beam is presented here. The peak enrgy 
is 5.19 kev. 
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The experiment was repeated with thicker crystals of 44 
lim and 320 Um thickness. Figures 16-19 are the plots of the 
data collection runs. The data in Figures 16 and 17 was 
taken with the crystal rotated to an angle of 0, =68.5°, 
while the data for Figures 18 and 19 was taken with the 
crystal rotated to 0, =104.8°. 

The target rotation angles used were obtained for 
maximum count rate and were determined by stepping the 44uUm 
portion of the crystal in 0.5 degree increments from 66.5° 
to 70° (<220> reflection plane) and from 101.6° to 104.8° in 
0.2 degree steps (<111> reflection planes). The highest 
count rates (photons per second) were found at 68.5° and 
104.8° respectivly. This would correspond to a Bragg angle, 
0, of 21.5° and 20.46° respectively. The count rate 
associated with each angle was 0.06 for <220> orientation 


and 0.06 <111> orientation. The angle of the _ Si(Li) 


detector remained fixed at 45°. 


D. PXR ANALYSIS 

To analyze the PXR data that had been collected, three 
independent methods were used. The first required the use 
the PCA software [Ref. 8}. The analysis included: 1) the 
determination of the gross area or total number of counts 
under the peak, 2) the number of background counts in the 


region of interest, 3) the gross peak height, 
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Figure 16. This plot was produced using the 44 wm 
Si <220> crystal orientation. The n=1 and n=2 
peaks were observed at 8.17 keV and 16.97 keV 
respectively. 
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Figure 17. This plot was produced using the 320 wu 
m Si <220> crystal orientation. The n=1 and n=2 
peaks were observed at 8.21 keV and 17.16 keV 
respectively. 
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Figure 18. This plot was produced using the 44 um 
Si <111> crystal orientation. The n=1,n=3, and 
n=4 peaks were observed at 4.56 keV, 15.31 keV, 
and 20.80 keV respectively. 
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Figure 19. This plot was produced using the 320 


um Si <111> crystal orientation. The n=l, 
and n=4 peaks were observed at 4.72 keV, 
keV, and 21.15 kev respectively. 
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n=3, 
15.62 


4) the peak FWHM, and 5) the centroid of the peak. For the 
software to provide this information, the region of interest 
(ROI) had to be defined. To define the region of interest, 
the peak was identified and then the leading and trailing 
edges were estimated. PCA software provided a fast way to 
return some initial values; however, the uncertainties in 
how the background was computed has left some question as to 
the accuracy of the results. 

The second method required the use of a spreadsheet 
[Ref. 9}. 

Before the data could be analyzed on the spreadsheet it 
had to be converted from binary to ASCII code. A program 
titled Convert written by Dr. George Rothbart was used. The 
program was provided by Dr. M. A. Piestrup of Adelphi 
Technologies Inc. With the data in ASCII format it was 
imported into the spreadsheet software. 

Using the spreadsheet, the areas--gross, net, and 
background--were computed. The observed resolution was also 
determined. As with the PCA software a region of interest 
was defined for each x-ray _ peak. The procedure for 
determining the ROI was different from that done using 
Reference 8. The background for the spectrum had to be 
determined. In the case of the <220> orientation spectras, 
the average background was a constant across the spectrun. 


For the <111> orientation a straight line was fit to 
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approximate the non-zero slope of the background. Once the 
background had been computed, it was subtracted from the 
spectrum. The ROI was defined as the area under the peak 
between the points where the tails went to zero. This ROI 
represented the net area. The gross area was then defined 
as the total area under the peak of the original spectrum 
between the points where the net area tails went to zero. 
The final method required plotting the spectrums and 
visually estimating the slope of the background with a 
ruler. The plots only allowed for estimating the resolution 
of the peak. The estimated values were within the 
“uncertainty of the spreadsheet values for 80% of the peaks. 
Tables 3 thru 8 show the comparison between the results 
generated by the PHA software and the spreadsheet software 
from Figures 16 thru 19. The resolution is not addressed in 
these tables since it will be addressed ina future tables. 
In Tables 3 thru 8 the following notation is used. E(n) 
is the x-ray peak energy of the nth peak. Ig (n) is the 
gross number of counts in the nth peak. Ip(n) represents 
the background counts in the nth peak. I(n) is the net 
counts in the nth peak. 
In another section of this chapter, the issue of 
intensity ratios will be addressed. When comparing the 
values of the intensities from Reference 9 and those from 


the PCA software, more than 80% of the intensity ratios were 
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within 15% of each other. Additionally, the PCA software 
provided no method for determining the uncertainty from its 
calculations. Consequently, the values of the intensities 
used were determined from the spreadsheet method. 

Another issue that will be addressed ina later section 
of this chapter is calculating the PXR linewidth. The 
measured linewidths of the observed spectra were used. The 
PCA software was not used as it gave linewidths that were as 
wide as the base of the peaks they represented. The 
uncertainties of the hand measurements are probably a more 
accurate estimate of the uncertainty than the uncertainties 
of the spreadsheet calculations. 

The uncertainty from using Reference 8, is not included 


in Tables 3 thru 8, as it was indeterminable. 


1. Intensity Analysis 

A transmission factor of less than unity reduces the 
X-ray intensity. The physical parameter to describe 
transmission is Lg, the absorption length in centimeters, 
which is dependent upon the absorption coefficient, UL. 
Table 9 is a list of wu and La for the theoretical PXR peak 
energies. The values of t in square centimeters per gram 
were determined using Reference 11. La, was computed using 


the following 


We = Se (18) 
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where p = 2.33g/cm3. 

The primary interest of this experiment was 
determining the effects on the higher order PXR at different 
crystal thicknesses. One way was to examine the relative 
intensities of the nth order peaks with respect to the first 
order peak. The intensity, I(n), (as discussed in Reference 


5) 1s 
2 
a 





I(n) « @p(n)Lal ze e Yel g, (@p(n))| (19) 
From equation (19) it is noted that the intensity is 
proportional to Lg- When Lg was greater than L, the crystal 
thickness along the particle beam path, with 
L = t/sinO,, (20) 
then the intensities of the higher order peaks were 
decreased. The converse is true. If L 1s greater than La, 
then the intensities of the higher order peaks will 
increase. Table 10 compares the values of L and Lg for the 
PXR data. 

The exponential term in equation (19) approaches 
unity as Lg becomes much greater than L. This accounts for 
the lack of higher order peaks in the 20 Um spectrums 
despite counting for two hours in the case of Figure 15 


(<111> orientation) and counting in excess of one hour for 


Figure 14 (<220> orientation). 
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Table 3. Comparison of PHA and spreadsheet (SS) analysis of 
data from 44 wm Si <220> orientaion. 


inst (PA) 2 (PHA) 1 (SS) 2 (SS) 
E(n) (keV) | 8.207+.035 |17.16+.04 |8.207+.035 | 17.16+.04 
E(n) /E(1) 2092-5013 2.094.013 
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Table 4. Comparison of PHA and spreadsheet (SS) analysis 
of data from 320 tim Si <220> orientation. 
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Table 5. Results of analysis by PHA of 44 um Si <111> 
orientation. 
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Table 6. Results of analysis by spreadsheet of 44 [im Si 
<111> orientation. 
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Table 7. Resuits of analysis by PHA of 320 Um Si <111> 
orientation. - 
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Table 8. Results of analysis by spreadsheet of 320 um Si 
<111> orientation. 
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Using equation (19), Table 11 was created to show how 


a change in x-ray peak energy can affect the intensity of 
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Table 9. Values of wt and La for theoretical parametric x- 
ray energy peaks. 
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Table 10. Comparison of L and La. 










the x-ray peak. The values listed our for a 20 Um crystal 


and a 44 jim crsytal. Table 11 reveals at low energies 
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cae! 2 can change as much as 5% while at higher 


energies it changes by only 0.1% 


mar [een [tm 











a> [p00 [as asz0.6 | ao9.20 
asx» [roe faa fatee Lno9.42 


(jim) 
Table 11. Example calculations of the effects’ on 
theabsorption factor in equation (19). 





Once the x-ray photons were generated, they passed 
through several media. The first was a one mil Kapton 
window (with a window diameter of one inch). Next the 
photons passed through a three centimeter air sank Finally, 
x-rays entered the Si(Li) detector via a two mil berrylium 
window. Tables 12 thru 14 list the attenuation coefficients 
for the PXR peak energies. Also found in these tables is 
the value of the absorption coefficient times the thickness 
of the media transitted, [px. Table 15 provides’ the 
attenuation factors of the Kapton window, the berrylium 
window and the air gap. 

Knowing the attenuation coefficient, the next step 
was to determine the intensity of the PXR prior to being 


transmitted throught the Kapton window. 
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px 

[<220>|e.a7 faa | 5. 784.08 | 020 5.0041 
<aso>|r6.97 | 773.004 |. 00274. 0005 
<220> |s.21 [320 | 5.708.08 |. 02024. 0042 
j<aso> jazae | |. 7534.004 | 0027+. 0005 
<ati> jase faa [33.848 | 2202. 024 
— 
eee 
320 












0035+. 0007 
.501+.002 .0018+.0004 
108.022 
<333>  |as.6e2_ | sf .9434.005 | .00334.0007 _ 
[<aaa> f2a.i5 | fe 6t.002 | .00174.0003 


Table 12. Attenuation coefficient for Kapton (C9905N9Hj0) 
with p = 1.42 g/cm3 and x=0.0025 + .0005cm. 
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8.17 1.064. 02 0096 
j<aso>|a6.97 | |. 263. 002 
1.05%. 02 
<aso> |a7ze | |. 260. 002 
<i> [4.se faa 5. 774.22 |. 05254.0063 _ 
5.32 | 299.002 | 0027 
<aaa> [20.80 | | sast.oor | 0002 
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2aee5 . 2 16eeooH .0002 


Table 13. Attenuation coefficient for berrylium. p = 
1.82g/cm3 and x = .0050cn. 





12Uncertainty was less than 2% except where noted. 


Using the following equation 
I(n) = I,(n)exp(-HPx) tot, (21) 
where 
(HPX) tot = 2(HPX) . (22) 
I, (n) is the intensity of the PXR after being created and 
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f<333>fas.aa | sation | -o06#.002 
<asa> [20.80 | | 72.003 | .0028+.0005 _| 
f<333>fas.62 | fa ast.o2 | 006¢.002 
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il ey SUB .685+.003 .0027+.0005 


Table 14. Attenuation coefficient for air. p = 0.001294 
X = 3.0+0.5cm. 





prior to exiting the scattering chamber. I(n) is the 
detector intensity. Rearranging the terms in equation (21) 
and solving for I, (n), equation (21) becomes 
Io(n) = I(n)exp(Upx) tot: (23) 
Table 16 is a compilation of the (Upx)tot values. Figure 20 
1s a plot of the attenuation factors of the Kapton window, 
the air gap, the beryllium window, and the product of the 
three attenuation factors versus x-ray peak energy. 
In Table 16, the final column lists the attenuation 


factors of the various peak energies. Not considered when 
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determining the factors was the efficiency of the Si(Li) 
detector. To determine the efficiency, the x-ray peak 
energy must be known. The peak energy determines whether 
the Be window or the detector thickness dominate to 
determine detector efficiency. For x-rays with peak 
energies below approximately 18 keV, only the beryllium 
attenuation factor is considered. For energies above 18 
keV, the detector thickness is the only contributing factor 
as the attenuation factor of the Be window is unity. [Ref. 
10) Figure 21 is reproduced from Reference 10, and is a 
plot of detector efficiency versus x-ray energy. 

The values in Table 16 for the total attenuation factor 
can be used to determine the detection efficiency except for 
the case of the <444> orientations. For the <444> 
orientations, Figure 10, must be included to account for the 
detector thickness contribution to the efficiency. grove 
these higher energy x-rays, the total attenuation factors in 
Table 16 would change from 0.995 to 0.945--approximately a 
5% difference. The corrected peak and area intensities are 
found in Tables 17 and 18. In the x-ray energy range 
examined, the detector has little bearing on the results; 
however, in ongoing experiments where PXR above 25 keV are 
being observed, the detector thickness can reduce the 


counting efficiency rapidly.[Ref. 10} 
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Br yd ern @ x 
(Kapton) (air gap) (berrylium) 


<z2o> faa fo.980fo.g6a sf o.990 
/fo.o97 Ss fos96 fog9a 
<aso> | i fo.997 fo.996 i fo.g98 
4g fo.ea7 i foeaa fonogo 
<333> | fo997 Ss fo.g94 fo.997 
/fong9e fo.997 fa 000 


Table 15. This table lists the attenuation factors for 
the x-rays going through the Kapton window, the air gap 
and the berrylium window. The differences in the <111> 
values of both crystal thicknesses is due to the energy 
differences between the two n=1 peaks. 





t (um) e| 


16.97 ___|.oo96+.0025 _|.990 
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Table 16. Total attenuation coefficient required for 
computing I,(n)- 
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Figure 20. Plot of attenuation factor for the 


Kapton window, the air gap, the Be window, and the 
product of the attenuation factors. 
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Figure 21. Calculated Full Energy Detection 
Efficiency Curve for Si(Li) X-Ray Detector. Loss 
in efficiency at high energies is due to detector 
thickness. At low energies only attenuation due 
to the beryllium window is considered. [Ref. 10] 
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prim) | rity [rym | tem /tya) 
239416 /255t18 


jaa | 239426 255418 
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Table 17. Comparison of the uncorrected peak intensities 
and peak intensity ratios with the corrected peak 
intensities and peak intensity ratios. 
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Table 18. Comparison of uncorrected area intensities and 
area intensity ratios with corrected area intenstites and 
area intensity ratios. 
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Tables 19 and 20 are comparisons of the corrected 
ratios with intensity ratios calculated from the theory. 
The theoretical ratios were reproduced from Reference 12. 
Note the discrepancies between the experimental peak ratios 
and the calculated ratios. The experimental area ratios 
compare more favorably with the calculated ratios. Figure 
22 is a plot of the ratios of the experimental peak ratios 
to the calculated intensity ratios versus peak number, n. 


Likewise, Figure 23 is a plot of the ratios of the 


— ex tose pean) ~~ wa 
| ep20> 
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Table 19. Comparison of the corrected peak intensity 
ratios with theory. Intensities are corrected for 
attenuation due the Kapton window, the air gap and the 
Be window. 


















experimental area ratios to the calculated intensity ratios 


versus peak number, Nn. 


oO] 


There was still discussion (Ref. 13) as to whether 
the calculated ratios were representative of the PXR peak or 
area intensity ratios. Further investigation is required to 


determine which case is correct. 


2. Resolution Analysis 
The secondary purpose of this experiment was to 
determine how the PXR resolution was affected by the change 


in PXR peak energy. 


I, (n)/I, (2) T,(m)/T, (1) (Exp area) 
(Calc) (exp area) / (calc) 





Table 20. Comparison of the corrected area intensity 
ratios with theory. Intensities are corrected for 
attenuation due to the Kapton window, the air gap, and the 
Be window. 

Recall from section B. of this chapter, that a 
straight line approximation was made to estimate the 


detector resolution (Fig. 8). Estimates of Iget were made 
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Figure 22. The plot shows a comparison of the 


experimental ratios divided by the calculated 
ratios for the peaks versus n, the peak order 
number. The numerical values are listed in Table 
19. 
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Figure 23. The plot shows a comparison of the 


experimental ratios divided by the calculated 
ratios for the areas versus n, the peak order 
number. The numerical values are listed in Table 


20. 
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and are listed in Table 21. Also in Table 21, are the 
values Tope of each PXR peak, and the derived value, Ipxr, 
of each x-ray peak. The latter values were calculated using 
equation (16) as modified by 

Téxr= Tope get (24) 
Table 22 is a comparison of the experimentally determined 
pyr and the Ipyp calculated from theory (Ref. 14]. Figures 
23 thru 26 are plots of the PXR linewidth versus PXR peak 
energy. The plots indicate that the linewidth increases as 
a function of energy. Another point of interest, the 
linewidths appear to be broaden with increasing crystal 
thickness for x-rays of the same energy. For each 
successive peak the value of AE/E decreases as E increases. 
AE is the linewidth of a given peak while E is the energy 
of the peak. Without more statistics, that is all that can 
be concluded about the linewidths. The issue as to how 
parametric x-radiation changes with energy requires further 
understanding. The starting point must be obtaining a 
better understanding of the detector and its resolution 


properties. 
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Table 21. This table is a list of the estimated detector 
linewidths and the PXR linewidths computed using Idet. 
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Table 22. Comparison of the experimentally determined PXR 
linewidths with the theoretically determined PXR 
linewidths from Reference 14. 
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Figure 24. PXR linewidth (keV) versus PXR peak 


energy for 44 Um <220> orientation. The slope m = 
0.0193 with an intercept b = .112 kev. 
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Figure 25. PXR linewidth versus PXR peak energy 


for the 320 Um <220> orientation. The slope is nm 
= 0.0156 and the intercept is b = 0.212 keV. 
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Figure 26. PXR linewidth versus PXR peak energy 


for the 44 wm <111> orientation. 


0.0242, the intercept is b = 0.1316 keV, 


correlation r = 0.9401. 
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The slope is m 


and the 
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Figure 27. PXR linewidth versus PXR peak energy 


for the 320 um <111> orientation. 


The slope is nm 


= 0.0188, the intercept is b = 0.320 keV, and the 


correlation is r = 0.9998. 
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IV. CONCLUSIONS AND RECOMMENDATIONS. 


The purpose of this experiment was to study parametric 
Xx-radiation from silicon crystals. PXR to order n=4 from 
Silicon were observed. Enhancement of the PXR from higher 
orders was determined to be due to the crystal thickness. 
The measurements required careful consideration of several 
experimental factors. 

First, an energy calibration was required to determine 
the values of the x-ray energy peaks. The calibration was 
also important in determining the resolution of the x-rays. 
The resolution of the PXR did not agree with the theory. 
However, the disagreement may be due to the assumptions made 
when determining the linewidths of the detector. Another 
explanation may be the method of calculating the theoretical 
linewidths. They were calculated without considering the 
detector apperature or the pipe that the x-ray pass through. 
These factors need to be reviewed to determine their affects 
on the x-ray linewidths. In future experiments, energy 
calibrations must be conducted prior to or following a data 
collection run. 

The second point was ensuring the angular position of 
the crystal was correct. The Bragg angle of 22.5° was the 
designated Bragg angle for the experiment. For the 20 Um 


crystal, this angle led to x-rays with peak energies that 
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concurred with theory. As for the 44 Um and 320 um 
crystals, the Bragg angle for the <220> orientation was 
21.5° while the Bragg angle for the <111> orientation was 
46. These angles produced bright x-rays at energies 
that were within ten percent of theory. Thus, despite the 
small angular deviations from the original Bragg angle, PXR 
was still observed and there is no other explanation for the 
x-rays. 

Finally, the intensity ratios were determined and 
compared with the theoretical intensity ratios. The peak 
intensity ratios deviated from the area intensity ratios by 
between 2% and 60%. The area ratio intensities concurred 
better with the calculated ratios than did the peak 
intensity ratios. Yet, the area ratios only agreed for the 
<220> orientation of the 44 Um and 320 um crystals. The 
<111> orientation differed by 12% to 24% for the 44 m 
cyrstal and from 35% to 70% for the 320 crystal. More data 
must be collected to determine if the correct intensity 
ratio to use is the peak ratio or the area ratio. From the 
results, preliminary indications suggest the area ratios are 
the correct ratios. The additional data is needed to also 
determine the effects on the intensity ratios when changing 
the Bragg angle by small increments. 

Having demonstrated the ability to produce low energy 


monochromatic x-rays, from a relatively small and low energy 
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machine, the next phase will be to work on methods to apply 
the x-rays to real world requirements. 

Before pursuing that endeavor more research is required 
into understanding how the resolution of the PXR is affected 
by the peak x-ray energy. Second, attention must be devoted 
to the understanding the electronics associated with the 
experiment and how the final results are influenced by such 
instruments. Finally, attention must be given to designing 
a more compatible target support ladder. On several 
occasions the ladder was a hindrance to the experiment 
rather than a support frame for the silicon crystal. With 
further research and investigation these problems will be 
solved. 

PXR is a new and promising field in x-ray generation. 
In time it will prove to be a reliable and safe way to 
produce bright x-rays in areas such as medical treatment 
facilities or small laboratories. PXR should be actively 


pursued. 
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